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2 ■ Abstract: 

Physical quantities in QCD are independent of renormalization scheme (RS), but that exact 
\ invariance is spoiled by truncations of the perturbation series. "Optimization" corresponds to 

^ . making the perturbative approximant, at any given order, locally invariant under small RS 

changes. A solution of the resulting optimization equations is presented. It allows an efficient 
algorithm for finding the optimized result. Example results for R^+e- = 3 ^ qf{l + TZ) to fourth 
order (NNNLO) are given that show nice convergence, even down to arbitrarily low energies. 
The Q = "freezing" behaviour, TZ = 0.3 it 0.3, found at third order is confirmed and made 
more precise; TZ = 0.2 ±0.1. Low-energy results in the MS scheme, by contrast, show the typical 
pathologies of a non-convergent asymptotic series. 



1 Introduction 



Renormalization, for physical quantities in massless QCD, amounts to eliminating the bare 
coupling constant in favour of a renormalized couplant. The precise definition of the renormalized 
couplant — the renormalization scheme (RS) — is in principle arbitrary, but at finite orders 
of perturbation theory the choice matters. It is well known that one can take advantage of 
this situation by allowing the renormalization scale, /i, to "run" with the experimental energy 
scale Q, but this familiar idea is vague and incomplete. What matters is not fi, in fact, but 
the ratio of /j, to A, itself a RS-dependent parameter — and, at higher orders, there are further 
sources of RS ambiguity. "Optimized perturbation theory" (OPT) [1] provides both a complete 
parametrization of RS dependence and the means to sensibly resolve these ambiguities by taking 
full advantage of Renormalization-Group (RG) invariance [2]. 

The aim of this paper is twofold: (i) to present a mathematical solution to the "optimization 
equations" of Ref. [T], allowing the optimized result to be found efficiently; and (ii) to show 
numerical results for R^+e- = <y[e^e~ — > hadrons)/(T(e+e~ — > iJ,~^fi~), updating the results of 
Refs. [31 mis] now that fourth-order calculations are available [6]. 

The RS-dependence problem remains a controversial topic. Our arguments have been set 
out in detail in Refs. [1] and [7]-|10j. Here we give only a brief exposition of the key idea 
behind OPT, the "principle of minimal sensitivity." This is the very general notion that, in 
any approximation method involving "extraneous" parameters (parameters that one knows the 
exact result must be independent of), the sensible strategy is to find where the approximant 
is minimally sensitive to small variations of those parameters. The unknown exact result is 
globally invariant, while the approximant is not; where the approximant has the right qualitative 
behaviour — local invariance (flatness as a function of the extraneous parameters) — is where 
one can have most confidence in its quantitative value. Many instructive examples testify to the 
basic soundness and power of this idea (see p].pT]-[T5]). 

In the present context, RG invariance tells us that physical quantities should be independent 
of fx and all the other "extraneous" parameters involved in the RS choice. That statement 
translates into an infinite set of equations (Eqs. (j2.16p below) that any physical quantity TZ 
must satisfy [T]. Perturbative approximations to TZ do not satisfy these equations, but we can 
find an "optimal RS" in which they are satisfied locally. We explain in this paper how these 
optimization equations can be solved efficiently. 

Calculations in QCD perturbation theory, in particular for /? and -Rg+e-; have progressed 
from leading order (LO) [16j, to next-to-leading order (NLO) [171118]. to next-to-next-to-leading 

^ The importance and generality of this idea was emphasized in [1] but earher authors had employed it in 
specific contexts, notably Caswell and Killingbeck [TT] in the anharmonic oscillator problem. 
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(NNLO or N^LO) ^ to now next-to-next-to-next-to-leading order (NNNLO or N^LO) 
\21\ |6]. (We will use the terminology "first order" for leading order (LO) and "(A; + 1)*^ order" 
for N'^LO.) Tribute should be paid at this point to the heroic efforts involved in these enormously 
complex calculations. The results allow us now to get a real sense of how QCD perturbation 
theory behaves, both in a fixed RS and with the optimization procedure. 

Perturbation series, in any fixed RS, are expected to be factorially divergent. However, it 
is possible that the optimized results converge, thanks to an "induced convergence" mechanism 
in which the optimized couplant shrinks from one order to the next [22] (see also \23\ [T2t [T5] ) . 
Our numerical results here are quite consistent with that idea. 

The plan of this paper is as follows: Section 2 reviews the mathematical consequences of RG 
invariance derived in Ref. [T]. Section 3 discusses finite-order approximants and explains how 
to "optimize" the RS choice based on the principle of minimal sensitivity. Section 4 solves the 
resulting optimization equations, in a general (fc + l)*'^ order, and Section 5 outlines an algorithm 
to determine the result efficiently. Section 6 applies this algorithm to obtain illustrative results 
for i?e+e- ™ second, third, and fourth orders and compares them with results in a fixed scheme, 
MS(/i = Q). The optimized results exhibit steady convergence — even down to zero energy, 
where the limit TZ = 0.3 ±0.3, found in third order [31 HI [5], is confirmed and made more precise; 
TZ = 0.2 lb 0.1. The fixed-scheme results show pathologies at low energies and we argue that 
these are a preview of what one can expect at higher orders at higher energies. Concluding 
remarks are in Section 7. Appendix A proves an identity mentioned in Section 4, and Appendix 
B discusses optimization in the fixed-point {Q — t- 0) limit, following Ref. |24j . 

2 RG invariance and its consequences 
2.1 RG equations 

We begin by reviewing the formalism introduced in Ref. [1]. First, we emphasize that OPT 
can only be applied to physical quantities (cross sections, decay rates, etc.). QCD involves 
many other objects - Green's functions, renormalized couplants, etc. - that are not physically 
measurable and are not RS invariant; for our purposes these are merely intermediate steps in 
calculating physical quantities and need not be discussed. In general, a perturbatively calculable 
physical quantity will have the form AiTZ+Aq with a leading-order coefficient Ai and, sometimes, 
a zeroth-order term Aq. The coefficients Aq and Ai, which carry dimensions of energy to the 
appropriate power, are RS invariant, so we may focus on dimensionless, normalized physical 
quantities TZ of the form: 

7^ = 0^(1 + ria + r2a2 + ...), (2.1) 
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where a = as/vr is the renormahzed couplant. The power p is usuahy 1 or 2 or 3, but need not 
be an integer. Generahy a physical quantity TZ is not a single quantity but rather a function of 
several experimentally defined parameters. We may always single out one such parameter with 
the dimensions of energy that we may call the "experimental energy scale" and denote by "Q." 
(It is needed only to explain which quantities are, or are not, Q dependent; the precise definition 
of Q in any specific case is left to the reader.) Note that Q must be defined in experimental 
terms and should not be confused with the renormalization scale /j,, which is defined only in 
terms of the technical details of the RS adopted. 

In the specific case of i?e+e- where, ignoring quark masses. 



^e+e- =3J]g,'(l+7^), (2.2) 



we have p = 1 and the e^e~ center-of-mass energy is a natural choice for Q. 

The fundamental notion of RG invariance f2] means that a physical quantity is independent 
of the renormalization scheme (RS). Expressed symbolically it states that 



^ dn dn 



da dlZ 



d{RS) d{RS) 

where the total derivative is separated into two pieces corresponding, respectively, to RS depen- 
dence from the rj series coefficients, and from the couplant a itself. A particular case of Eq. 
()2.3p is the familiar equation expressing the renormalization-scale independence of TZ: 

d 



+ I3{a)^]n = 0, (2.4) 



where 



/3(a) = = -6a^(l + ca + C2a^ + csa^ + ...). (2.5) 

dpi, 

The first two coefficients of the /3 function are RS invariant [25] and, in QCD with massless 
fiavours, are given by [T7]: 

(33 - 2n/) _ 153 - 19n/ 



6 ' 2(33 -2n/)' 

When integrated, the /3-function equation can be written as: 



(2.6) 



da! ^ i'^ du! , , , 

+ C= / -^i = ln(/i/A). (2.7) 



where C is a suitably infinite constant and A is a constant with dimensions of mass. The 
particular definition of A that we use corresponds to choosing [1] 

C= / — 97 r 2.8 
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(where it is to be understood that the integrands on the left of ()2.7|) are to be combined before 
the bottom hmit is taken) . This A parameter is related to the traditional definition [26\ [27] by 
an RS-invariant, but n j-dependent factor [Tl [28]: 

ln(A/A) = (c/6)ln(2c/6). (2.9) 

The A parameter is scheme dependent, but the A-parameters of two schemes are related by 
the Celmaster-Gonsalves relation [29]. If two schemes are related, for the same value of /j,, by 



a' = a{l + via + ...) (2.10) 

then 

ln(A7A) = vi/b. (2.11) 

This relationship is exact and does not involve the V2,V3,... coefficients. Thus, the A's of 
different schemes can be related exactly by a 1-loop calculation. Hence, the A parameter of any 
convenient "reference RS" can be adopted, without prejudice, as the one free parameter of the 
theory, taking over the role of the "bare coupling constant" in the original Lagrangian. 

From Eq. (j2.7p it is clear that a depends on RS only through the variables fi/A and C2, C3, . . ., 
the scheme-dependent /3-function coefficients. The coefficients of TZ can depend on RS only 
through these same variables — the RG-invariance equation, (j2.3p . could not be satisfied oth- 
erwise. Therefore, these variables provide a complete RS parametrization, as far as physical 
quantities are concerned [T]. Thus, we may write: 

o = a(RS) = a(r, C2, C3, . . .), (2-12) 

where 

T = bln{fj,/A). (2.13) 

The T variable is convenient and helps to emphasize the important point that /i itself is not 
meaningful because of the scheme ambiguity represented by Eq. (|2.10p ; only the ratio of to A 
matters. 

The dependence of a on the set of RS parameters r and cj is given by [l]: 

fin 

- = /3(a)/5, (2.14) 



13 j {a) = -hl3{a) / dx (2.15) 



da x^^"^ 
— = /3j{a) = -b/3{a) / dx —— 
dcj Jo P{x) 

where the first equation is just the /3-function equation in new notation, and the second follows 

by taking the partial derivative of Eq. (j2.7p with respect to Cj, holding r and the other {i 7^ j) 

constant. 
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The symbolic RG-invariance equation ()2.3p can now be written out explicitly as the following 
set of equations [Ij: 

(2.16) 

d_ 

da 



dn 


(- 






an 


(- 


dcj ~ 


ydcj 



+ /3,(a)— ) 7^ = j = 2,3,..., 



2.2 The /3j{a) and Bj{a) functions 

The functions, defined in Eq. (j2.15p . begin at order a^~^^ so it is convenient to define Bj{a) 
functions whose series expansions begin 1 + 0(a): 

B,{a)^^-^^,{a). (2.17) 

For j = 1 it is natural to define 



oo 



Bi{a) = B{a) = 3^ = 1 + ca + caa^ + . . . = ^ c,a\ (2.18) 
" i=0 

with the convention that cq = 1 and ci = c. Equation ()2.15p can then be re-written as 

B,{a) = ^^-^B{a)I,{a), (2.19) 



where 



(Note that this formula for Bj(a) even holds for j = 1 if the r.h.s. is interpreted as the limit 
J — )• 1 from above.) 

The Bj{a) functions have power-series expansions whose coefficients we write as 

00 

= (2.21) 

with Wq = 1. The other W- coefficients are fixed in terms of the Cj's [Ij. Differentiating 
Eq. (|2.15p . or equivalently by requiring commutation of the second derivatives, d'^a/drdcj = 
d'^a/dcjdr, leads to [T] 

I3'j{a)/3{a) - /3'{a)/3j{a) = -ba^+^, (2.22) 

where here the prime indicates differentiation with respect to a, regarding the coefficients cj 
as fixed. From this differential equation it is straightforward to show that the W-^s satisfy the 
relation 

i 

(f + j - 1 - 2m)cmWi^ = (j - l)dio (2.23) 

m=0 



5 



for i = 0, 1, 2, . . . and j = 1,2, . . .. In the special case j = 1 one has = Ci and the above 
equation reduces to 

i 

^ (i - 2m)cmCi-m = 0, (2.24) 



m=0 

which is true identicahy, since the left-hand side is 



^ ^ {i ^)CmCi—m ^ ^ nT'CmCi—m (2.25) 

and the first sum, by changing the summation variable from m to n = i — m, is seen to cancel 
the second. 

2.3 The pn invariants (p = 1 case) 

The RG equations (|2.16p determine how the coefficients rj of TZ must depend on the RS variables 
{r, Cj}. To show explicitly how this works we specialize to the p = 1 case, where 

-JZ = a{l + ria + r2a^ + . . .) (2.26) 

and write out the lowest-order terms to obtain 



a^-^ + a''-^ + - a^(l + ca + ...){! + 2ria + ...) = 0, (2.27) 
a^l^ + a^^ + ..A +a\l + W^a + . . .)(1 + 2na + ...) = 0, (2.28) 

OC2 OC2 J 

and so on. (In fact, the coefficient is zero.) Equating powers of a, one sees that ri depends 
on r only, while r2 depends on r and C2 only, etc., with 

dri 

^ = 1, (2,29) 

If—' 

etc.. Upon integration one will obtain as a function of r, C2, . . . , Q plus a constant of integration 
that is RS invariant. Thus, certain combinations of series coefficients and RS parameters are 
RS invariant [1]. The first two are 

PiiQ)=T-n, (2.31) 

P2 = C2 + r2 - cri — rj. (2.32) 

The first invariant, pi{Q), is unique in being dependent on the experimental energy scale, 
Q. A calculation of the coefficient ri, in some arbitrary RS, yields a result of the form 

ri = 61n(/i/Q) + ri,o, (2.33) 
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whose ^ dependence indeed conforms with Eq. ()2.29p . For dimensional reasons, the /i and Q 
dependences are tied together; the ri calculation does not "know" what boundary condition 
will later be applied to the /3-function equation, so the parameter A cannot explicitly appear. 
Similarly, the higher coefficients r2, . . . depend on ln{fj,/Q), but not on or Q separately. Hence, 
for the invariants p2, ps, ■ ■ ■ the cancellation of p dependence also implies the cancellation of Q 
dependence. However, pi{Q) is different because its definition explicitly involves r, and we find 

Pi(Q) = b\n{p/K) - {bHp/Q) + ri,o) 
= 61n(Q/A)-ri,o 

= 61n(Q/A7^), (2.34) 

where A is the A-parameter in the RS in which the ri calculation was done, and K-ji is a 
characteristic scale specific to the particular physical quantity IZ. We can regard the last step as 
a Celmaster-Gonsalves relation, (j2.1ip . that relates At^ back to the theory's one free parameter, 
the A of some reference RS. 

Some convention must be adopted to uniquely define the higher invariants pj (for j > 2) 
because, of course, any sum of invariants is also an invariant. For example, one might quite 
naturally add some multiple of to Eq. (j2.32p . Indeed, the tildes over the /jj's are included to 
distinguish them from an earlier definition One convenient definition is as follows |31] . For 
any given physical quantity TZ one can always define a RS (known either as the "fastest apparent 
convergence" (FAC) or "effective charge" [30] scheme) such that all the series coefficients 
vanish in that scheme, so that TZ = apAc(l + + + . . .). Since the /3 functions of any two RS's 
are related by 

da' da' da da' ^, , , , 

/3'(a') ^p-r = -r^'T- = :r^(«)' (^.35) 

dp da dp da 

we must have 

/^FAc(7^) = ^/3(a). (2.36) 
da 

The pn invariants can be defined to coincide with the coefficients of the FAC-scheme /3 function: 

n=0 

where po = 1 and pi = c (not to be confused with the independent invariant pi{Q) = r — ri). 
Re-arranging this equation as 

B{a) = Y^p^a-(-) ^ (2.38) 

n=0 ^ ^ da 
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and equating powers of a we obtain 



Cj — Pi Cj-i 



i=0 



da 



(2.39) 



where C„[-F(a)] means "the coefficient of a" in the series expansion of -F(a)." 
The first few invariants are hsted below (for p = 1): 



Pi 
h 
h 



c, and 
C2+r2- cri 



Pi{Q) =T-ri, 



1' 



(2.40) 



C3 + 2r3 — 2c2ri — 6r2ri + erf + Arf. 



(Note that our earher papers used a different convention, with = P2 — \c and p'^ 

2.4 The pn invariants (general p) 

For general P the first few invariants are 



.old 



PS- 



pi = c, and 

r2 cri 



Pi{Q) = '^-7' 



P2 



P3 



C2 + 



P + 1 

2p2 



IJ 



(2.41) 



, 2r3 2c2ri „ /P + 2 
C3 H 2 



r2ri + 



err 



+ 



2(p + l)(p + 2) 



p p V y p^ V 3p3 

The generalization of Eq. (j2.39p to general P is given below in Eq. (|5.1|) . An explicit inverse 
formula giving the p's in terms of the Cj and rj coefficients can be found [31] : 

n da) 



1 (-R}/^^ 



(2.42) 



3 Finite orders and "optimization" 
3.1 Finite-order approximants 

So far our discussion has been formal and the results have been mathematical theorems. Now we 
need to discuss finite-order approximants. At this point matters inevitably become controversial, 
because any approximation (unless it uses rigorously proven inequalities that bound the exact 
result) is necessarily a gamble; one is trying to guess at the exact result based on incomplete 
information. The issue is how best to use all available information. 

The first point to make is that two truncations are involved, for TZ and for (3. (We need (3 
in order to relate a to the A parameter of some reference scheme.) Thus, the {k + l)*'^ order, or 
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(next-to)'^-leading order (N'^LO) approximant is naturally defined with both TZ and /3 truncated 
after k + 1 terms: 

TeC^+i) =aP(l+ria + ...+rfca'^), (3.1) 
where a here is shorthand for a^'^"*"^^ the solution to the int-/3 equation with /3 replaced by 

= _ba'^(^i + ca + . . . + Cfco'^). (3.2) 

It is straightforward to check that the order of the error term TZ — 7^^^+^) is determined by 
whichever truncation, of TZ or /3, is the more severe, so it is natural to use the same number of 
terms in each [U [7] . 

In a fixed RS (with the RS choice also entailing a choice of fi), the first step will be to find 
the value of a in that RS by solving the integrated /3-function equation, Eq. ()2.7p . That equation 
can be re-written (with r = 61n(/i/A)) in the form: 

T = K{a) = K^^\a) - A{a), (3.3) 

where 



x^[l + cx) a 



ca 



1 + ca 



(3.4) 



and 

In {k + l)*'^ order B{x) is replaced by B^'^^^^x) = 1 + ex + . . . + CkX^. Hence A(a) would vanish 
in second order, so that K^'^\a) is indeed the second-order approximation to K{a). 

We shall refer to Eq. (j3.3p as the "integrated /3-function equation" or "int-/3 equation." It 
should be solved numerically — to an accuracy comfortably better than the expected error in 
the final result (see discussion in Section [6]). To use an analytic approximation, such as a trun- 
cated expansion in inverse powers of ln(;U^/A^) [261 127j. would introduce another uncontrolled 
approximation and create another source of ambiguity |321 [8] , namely dependence on how pre- 
cisely the A-parameter is defined (e.g., the choice between A and the more conventional A; see 
Eq. (j2.9p ). This is a wholly avoidable ambiguity and it is sensible to avoid it. 

3.2 Optimization in low orders (p = 1 case) 

While the exact TZ is RG-invariant, the finite-order approximants are not, since the truncations 
spoil the cancellations in the RG equations (j2.16p . If TZ in those equations is replaced by 
is not zero but is some remainder term 0{ciF^^^^). As explained in the 
introduction, the idea of "optimized perturbation theory" is to choose an "optimal" RS in 
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which the approximant T^^'^^^) is locally stationary with respect to RS variations; i.e., the RS 
in which 7^^^+^) satisfies the RG equations, ()2.16p . with no remainder: 



O-jZik+i) 



y dr 



+ 



/3(a) 97^('=+l)' 



da 



opt. RS 



+ /3,(a) 



da I 

/ opt. RS 



0. 



0. 



(3.6) 



We assume here that the QFT calculations of the TZ and /3-function coefficients up to and 
including and have been done in some (calculationally convenient) RS. From those results 
we can compute the values of the invariants Pi{Q) and pi = c and p2, ■ ■ ■ ,Pk- Our optimized 
result will be expressed solely in terms of those invariants. 

To see how this works let us consider the second-order (NLO) approximant. (For simplicity 
we set p = 1 for the remainder of this subsection.) 



7^(2) = a{l+ria), 



(3.7) 
(3.8) 



where a here is short for a^'^\ the solution to the int-/3 equation (j3.3p with /3 replaced by /S^^^ 

1 



K(2)( 



a ) = — h c In 
a 



ca 



1 + ca 



(3.9) 



Since TZ^'^^ depends on RS only through the variable r, only the "j = 1" equation in Eq. 
above is non-trivial. Thus, the optimized TZ^"^^ is determined by a single optimization equation: 



^ - a2(l + ca)(l + 2fia) = 0. 



(3.10) 



(Overbars are used to indicate the value in the optimum RS.) As discussed in subsection 12.31 
the a^ terms must cancel in any RS, which fixes ^ = 1, leaving 



l-(l + ca)(l + 2fia) 



0. 



(3.11) 



This determines the optimized coefficient ri in terms of the invariant c and the optimized 
couplant a: 

1 c , , 

ri = (3.12) 



2 1 + ca 

But fi is related to f by the definition of the Pi(Q) invariant, Eq. (12.310 : 



PiiQ) ^T-n. 



(3.13) 
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Eliminating ri between these last two equations and substituting into the second-order int-/3 
equation, (j3.9p . gives 



, , ca \ 1 f ca 
1 + CO In + 



l + ca/ 2 \1 + ca 



PiiQ)- (3.14) 



If the values of the invariants c and Pi{Q) are known, as we assume, then we may numerically 
solve this last equation to obtain a. Substituting back in (I3.12p we can find fi and hence obtain 
the optimized approximant TV-"^^ = a(l + fio). 

Note that the only approximations made here are the truncations of the TZ and (3 series 
in Eqs. (j3.71 13. Sp . which define the second-order approximant is some general RS. We do not, 
for instance, approximate Eq. (j3.12p as fi ~ — (which corresponds to the PWMR [33] 
approximation, discussed later). Also, we will need to solve Eq. (|3.14p numerically 

We now turn to third order. The third order approximant is defined by 

/3(3) = -60^(1 + ca + C2a2), (3.16) 

where now a is short for a^^\ the solution to the int-/3 equation with /3 replaced by /S^^). Since 
7^('^^ depends on RS through two parameters, r and C2, there are two optimization equations 
coming from Eq. (|3.6p . These correspond to suitably truncated versions of Eqs. (j2.27l I2.28p . 
Using Eqs. ^CTM [230]) they reduce to: 

1 + (2fi + c) a - B^^\a) (l + 2fia + 2,7=2(1^) = 0, (3.17) 

1 - Bf\a) (1 + 2fia + Sfsa^) = 0. (3.18) 

Here B^^\a) = {l + ca + C2a^) is the B{a) function truncated at third order in the optimum 
scheme. The other function B2 (a) is not a polynomial; it is obtained from Eqs. (j2.191 12.20p 
with B{x) replaced by (1 + cx + C2X^). 

These optimization equations, together with the definitions of the invariants p2 and pi{Q) 
and the int-/3 equation, fully determine the optimized result, TZ^^^ [HIS]. The new algorithm of 
Section [5] provides a more efficient route to the result than the method in Ref. [5] . 

Note that the "optimum RS" evolves from one order to the next; for instance fi at third 

(2) 

order is not the same as fi at second order (so, strictly we should have distinguished fl and 
f^^^ in the above). 
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3.3 The optimization equations 

We now write down the optimization equations at some general, {k + 1)*'', order [l]. To treat a 
general power p, it is convenient to define 

S = ^^, (3.19) 

pa^~i da 



whose series expansion 



has coefficients 



S = l + sia + S2a^ + ... (3.20) 



m + P 



, ... (3-21) 
p , 

Using the Sm coefficients absorbs all the P dependence in the analysis of the next section. 
(However, p will re-appear later when we need to combine those results with the P-dependent 
Pn invariants.) 

As we saw in subsection 12.31 (in the p = 1 case) all terms in the RG equations up to and 
including 0{aF~^^) must cancel automatically in any RS. In the "j = 1" optimization equation of 



Eq. m, the 



term is a polynomial which must cancel the first k terms of i^Mi^Eii^l^ 



b da 

A similar observation applies to the other optimization equations of Eq. (j3.6p . Hence, we may 
reduce the optimization conditions to 

5]'+^) (a)5(^+i) (a) - T,_,- [^f {a)S^''+'^ (a)] = 0, (3.22) 



for j = 1,2, ... ,k, where the notation T„[F(a)] means "truncate the series for F{a) = Fq + 
Fia + . . . immediately after the a" term" (i.e., T„[F(a)] = Fq + Fia + . . . + Fna"-.). 

Note that the B^''^^\a) functions are defined by Eqs. (l239l [2:20]) with B{a) replaced by 
the polynomial B^'^'^^^a) = 1 + ca + . . . + Cka^. 



4 Solution for the optimized coefficients 

4.1 Definition of the Hi(a) functions 

In this section it is implicit that all quantities are in the optimum RS at (k + l)^^ order; overbars 
and (fc+i) superscripts will be omitted. Also, we make the convention that 

ro = So = Co = 1, and ci = c. (4.1) 

Next ~ for reasons that will become clear in the next subsection - we define some functions 
Hi(a), . . . , Hjt{a) that are combinations of the Bi{a), . . . , B^^a) functions: 

g,(a)^^c,a^ j ^i3,+,(a) i = (1), 2, . . . , fe. (4.2) 
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For i = 1 this definition, as it stands, is ambiguous; it should be interpreted as 

Hi{a) = Bi{a) - cja^Bj+iia), (4.3) 

corresponding to 

It is also convenient and natural to define 

Hoia) = 1 and Hk+i{a) = 0. (4.5) 

The i/'s are defined as combinations of the -B's. It turns out that there is a simple formula 
for the inverse relationship, giving the i?'s as combinations of the H^s. 
Lemma: 

k-j 

Bj{a) = J2 W^amj+g{a) j = 1, . . . , A;, (4.6) 

where the W- coefficients are those of the series expansion of Bj{a), Eq. (j2.2ip . (One might 
describe this result as follows: Take the power series for Bj{a) and truncate it after the a''~^ 
term. Now re- weight each term, replacing a*^ by a'^Hjj^q{a), and the result is the full series for 



Proof: We first treat the cases with j ^ 1. Using the definition of the -ff's, Eq. (j4.2p . the 
r.h.s becomes 

Reorganizing the double sum by defining n = q+ p converts this expression to 

E ^^T^ + 2v)o,WU. (4.8) 

n=0 "^J^^ ^1 p=0 

The inner sum reduces to (j — l)(5no by virtue of Eq. ()2.23p . Thus, only the n = term of the 
outer sum survives, the (j — 1) factors cancel, and one is left with just Bj(a), as claimed. 
In the case j = I the result to be proved, Eq. (j4.6p . becomes 

fc-i 

Bi{a) = Y,Cga''Hg+i{a). (4.9) 
Using Eq ()4.3p for Hi{a) and Eq (14. 2p for the other H's, the r.h.s. becomes 

k—l k—l k—q—l , 

j = l q=l j=0 W+J^ 
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Reorganizing the double sum by defining n = q + j yields 

a: 



Bi{a) - Cja^Bj+i{a) + ^ — fi„+i(a) ^(2g - n)cgCn-q (4.11) 

j = l n=l " \q=l J 

The inner sum, in parentheses, after adding and subtracting a q = term becomes 



'Y{2q - n)cgCn-q + ncn = ncn (4.12) 

9=0 

since the full sum vanishes, as noted in Eq. (|2.24p . Thus, the two series terms in (j4.1ip cancel 
leaving just Bi{a), as claimed. 

4.2 Formula for the optimized coefficients 

We are now ready to state the main new result; an exact, analytic expression for the optimized 
Sm (and hence the r^) coefficients, for m = 0, 1, . . . , /c, in terms of the (optimized values of) a 
and the /3-function coefficients C2 , . . . , : 
Theorem 

The optimization equations ()3.22p are satisfied by 

Sm = ^ , . (Hk-mia) - gfc-m+i(a)) , m = 0, 1, . . . , fc. (4.13) 

Bk{a) 

Proof: hi the case j = k, where the To[. . .] term in (j3.22p is just unity, the optimization 
equation reduces to 

We first prove that this equation is satisfied. Substituting Eq. (I4.13P into the series for S gives 
S=Y^ Sma"^ = Yl i^k-M - Hk-m+i{a)) . (4.15) 

m=0 ' m=0 

The i?'s cancel in pairs leaving 

S = {H, - Hu+i) = (4.16) 

since we defined Hq = 1 and Hf^^i = above. 

Using this result and writing out the truncated-series term explicitly, the remaining opti- 
mization equations can be re-written as 

f7^ = i:«^E--^t™ J = l,...,k-l. (4.17) 
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We now need to prove that these equations are satisfied by Eq. (I4.13p . The r.h.s. becomes 

a"™ 



E E ^'-m^ {Hk-m{a) - Hk-m+i{a)) . (4.18) 

i=0 m=0 



Reorganizing the double summation, defining q = i — m and thereby replacing i with m + q 
yields 

k-j k-j-q 

TfyTE^'^' E {Hu~M - Hk-m+i{a)) . (4.19) 

g=0 m=0 

The inner summation reduces to Hj^q{a) since the iJ's again cancel in pairs (and -fffc+i = 0)- 
Thus, the r.h.s of (|4.17p reduces to 

(4.20) 

where the last step uses the Lemma, Eq. ()4.6p . and produces the l.h.s. of ()4.17p . completing the 
proof. 

4.3 An identity and the PWMR approximation 

It is worth noting the following "complete sum" identity (proved in Appendix B) 

E^i"' ( ■ + jl| )^^+^-(°) = ^ i = (l),2,...,fc, (4.21) 

with the i = 1 case interpreted using (j4.4p . This identity reveals a remarkable property of the 
//i(a)'s, which are defined as a "partial sum" (over j = 0, . . . ,k — i) of the same terms. Hence, 
we can write 

H,{a) = l- 9«^'( '~'^.l! )^-^+i(«) i = {l),2,...,k, (4.22) 

which, unlike the Hi definition, involves Bj's with j greater than k. Since the Bj^s all start 
1 + 0{a) we see that the series for {Hi{a) — 1) begins only at order a^^*+^: 

Hi{a) - 1 = tl^^Ck-^+la''~'^' (1 + 0(a)) . (4.23) 

Substituting this result into Eq. ()4.13p quickly leads to the result 

Sm = — ^ — c^ + 0(a). (4.24) 

This result was first obtained — in a quite different manner — by Pennington, Wrigley, and 
Minaco and Roditi (PWMR) [33]. The resulting PWMR approximation can be useful when 
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a <^ 1. (In fact, one also needs a <^ a* if a "fixed point" a* exists; see Appendix B.) Inserting tlie 
above equation into the definitions of the pn invariants, one can find the PWMR-approximation 
Sm's in terms of the invariants. For k = 2 one finds 

3 

si~0, S2~--P2, (4.25) 

while for A: = 3 one finds 

si^lc, S2«-^P2-^c^ ss^-2p,-lp2C-^c\ (4.26) 

These results provide a useful starting point for the optimization algorithm described in the 
next section. 



5 Optimization algorithm 

The optimization problem at {k + 1)*'^ order involves 2A; + 1 variables, namely, a, r, C2, . . . ,Cfc, 
and ri, . . . , Tfc. These are connected by 2/c + 1 equations, namely, the int-/3 equation, the k opti- 
mization equations, and the k formulas for the invariants Pi{Q) and p2, ■ ■ ■ , Pk, whose numerical 
values we assume are given — in the case of Pi{Q) the numerical value will depend on the value 
of Q being considered. We shall use a,C2, ■ ■ ■ ,Ck as the principal variables. The solution to the 
optimization equations then explicitly determines the coefficients ri , . . . , in terms of these 
principal variables. The int-/3 equation explicitly fixes r in terms of the principal variables. The 
PiiQ) = T — ri equation can be used at the end to relate a to Q, so the remaining task is 
to use the formulas for p2, . . . ,pk to determine, self-consistently, by some convergent iterative 
procedure, the C2, . . . , Cfc variables. One such algorithm is the following (recall Sm = ( "^+^ ) r^): 

(1) Choose a numerical value for a. 

(2) Make an initial guess for numerical values of si, . . . ,Sk 
(e.g. use the PWMR approximation). 

(3) Find cj values from the invariants using Eq. (j2.39p generalized to any P: 

j 

Cj = Pi Cj_j 

i=0 

where C„[-F(a)] means "the coefficient of a" in the series expansion of F{a)." 

(4) Hence construct the B{x) function as X]j=o'^i^"' ^^'^ then obtain B(a) by sub- 
stituting X = a and the Bj{a)^s {j = 2,...,k) by numerical integration of their 
definition: 

B,ia)=^-l-^Bia)£dx^. (5.2) 



V / nda 
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(5) Now find the /7j(a)'s from their definition: 



j=o ^ ^ 



(5.3) 



and hence obtain new values for the Sm coefficients from 



a 



—m 



Bk{a) 



{Hk~m{a) — Hk-m+i{a)) ■ 



(5.4) 



(6) Iterate from step 3 until the results converge to the desired precision. 



(7) Finally, use 




(5.5) 



(from the int-/3 equation r 



= K^^~^^\a), combined with pi{Q) 



T — Ti/p) to find 



the value of Q that corresponds to the chosen a value. 

(8) One can then repeat the whole procedure with different initial a values to cover 
the desired range of Q values - or to home in on one particular Q value, if desired. 

Various details of this algorithm can be refined. In particular one can try to avoid numerical 
problems with the large cancellation between the two H^s in ()5.4p . by constructing Hi{a) — 1. 
This algorithm appears to be quite robust and efficient provided that B{a) is not too small; 
i.e., provided we are not close to a fixed point. In that case one can find one of the iterations 
generating a negative B{a) which makes the integrals in (j5.2|) blow up. A simple cure for this 
problem is to is to use a modified algorithm that works with a fixed B{a) rather than a fixed a. 
That is, one sets a fixed value of B{a) = Bq at step 1 and then in step 4, having constructed 
the new B{x) function, one solves for a new a from B{a) = Bq. Also, near to the fixed point the 
PWMR approximation is not a good starting point and a better one is to use the fixed-point 
result instead (see Appendix B). 

6 Numerical Examples for Re+e~ 
6.1 Procedure 

We turn next to illustrative numerical results for a specific physical quantity, namely 



o"tot(e''~e — >■ hadrons) 



(6.1) 



cr(e+e —?■ i^i^ fj, ) 
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at a centre-of-mass energy Q. We shall neglect quark masses and consider nj flavours of quarks 
with electric charges (i = 1, . . . ,?^/)• If there were no QCD interactions then Rf,+e- would 
equal the parton-model result, 3^g?. Including perturbative QCD corrections we have 

^e+e- =(3 5^(?,')(l+7^), (6.2) 

where 7^ is a normalized physical quantity whose perturbation series has the form 

n = a{l + ria + r2a^ + ...). (6.3) 

Previous discussions of this quantity to third order in OPT O HI [S] can now be extended to 
fourth order thanks to recent results from Baikov aZ p]. Most importantly (see Example 7 
below), we find that the infrared fixed-point ("freezing") behaviour is confirmed and made more 
precise in fourth order. There is, though, no need to update the phenomenological conclusions 
of Ref. [5] because the changes are well within the uncertainties discussed there. Therefore, we 
shall ignore various phenomenological issues (the effect of quark masses, the matching of A's at 
flavour thresholds |34], Poggio-Quinn- Weinberg smearing |35j, experimental uncertainties, etc.) 
that were discussed in Ref. [5j. 

Our focus here will be on the apparent convergence, or otherwise, of results from second to 
third to fourth order. We shall compare the optimized results with those in a fixed scheme, the 
"modified minimal subtraction" (MS) scheme with the renormalization scale /i chosen equal to 
the center-of-mass energy Q. (Properly speaking, we should denote this scheme as MS(/i = Q), 
but for brevity we leave the fi = Q specification understood.) We will essentially presume that 
the value of A^g is known from fitting other experimental data. However, to avoid committing 
to any specific value, we label our examples, not by Q, but by the ratio of Q to A^. At each 
order we proceed as if only the coefficients to that order had been calculated. 

To obtain the MS results, the first step is to evaluate the numerical value of the r parameter 
of the MS scheme: 

r^ = 61n(Q/AMs). (6.4) 

Then, one must numerically solve the int-/3 equation, p.3p . At second order, where B{a) is 
approximated by 1 + ca, this equation is 



= i^(2)(a) = i + cln 
a 



ca 



1 + ca 



(6.5) 



(2) 

With the resulting a, one then evaluates = a{l+r^^a). At third order one must numerically 
solve 

r^ = Kg(a) = i^(2)(a)-Ag(a), (6.6) 
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where 

Jq \B{x) l + cx J 
with, in this case, B{x) approximated by 1 + cx + C2x'^ with C2 = c^^. (A convenient approach 
is to first evaluate A at some initial value of a; then numerically solve Eq. (j6.6p . with A treated 
as a constant, to obtain a new a; and then iterate.) With the resulting a one then evaluates 
= a(l + r^^^a + r^^a^). At fourth order the procedure is the same, except that one now 
includes a c^^ term in B{x) and an r^^ term in TZ. 

To obtain the optimized results one first needs to calculate the numerical values of the 
invariants. At a given Q/A^g one can find from Eq. (j6.4p and then obtain pi{Q) as r"^^ — rj^^. 
Numerical values for the p2 and ps invariants (which are Q-independent) are similarly obtained 
by evaluating their definitions in Eq. (|2.4U|) using the MS-scheme and Cj coefficients. The 
optimized result to second order is obtained from Eqs. (j3.14p and (|3.12p . At higher orders one 
uses the algorithm described in the preceding section. 

At each order one wants, not only a result for IZ but also an estimate for its likely error. There 
is no rigorous way of doing this. However, it is reasonable to expect the "apparent convergence" 
of the series (i.e., the behaviour of the terms that have been calculated) is some sort of guide. We 
shall adopt the common practice when dealing with asymptotic series of viewing the magnitude 
of the last calculated term, | r^a'^"'"^ |, as the error estimate. We do this both for the MS and the 
optimized results. The change in the TZ results from one order to the next is another indicator; 
it seems quite consistent with our error estimate. 

We will give two sets of examples; one set at moderately high energies, and the other at low 
energies. For the first set of examples the appropriate number of flavours is = 5 {u,d,s,c,b 
quarks), while for the second set it is rij = 2 {u,d quarks only). 

6.2 High-energy examples 

For Uf = 5 the /3 function's leading, RS-invariant coefficients are 

, 23 29 , , 

6=—, c=— . 6.8) 
6 23 ^ ^ 

In the MS scheme its next two coefficients are [19^ [2T] 

9769 

= ^ = 1.474789, (6.9) 
26017 11027 

^ 31104 1242 ^ ^ ^ 

where ("(s) is the Riemann zeta function. The MS coefficients in TZ{e~^e~) are \18 \ \20 \ [6] 

rp = 1.409230, rp = -12.80463, rp = -80.43373. (6.11) 
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(The exact values, involving vr^, C(3), C(5) and [6] were used in our calculations.) Inserting 
these values in Eq. ()2.40p yields 



P2 = -15.0926, p3 = -33.2216. (6.12) 

Our first two examples update those in Table II of Ref. [5J. Results for {k + I)*'* order 
(k = 1,2,3) in both the MS and optimized schemes are presented in the tables and figures 
below. 

At these energies the perturbation series seems well behaved. The MS results are quite 
satisfactory but the optimized results offer greater precision, with smaller expected errors that 
tend to shrink more rapidly with increasing k. It is noteworthy that while Oms slightly increases 
with k, the optimized couplant a shrinks, consistent with the "induced convergence" scenario of 
Ref. i22j. 
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Example 1: Q/Ams = 340. 



Order 




T^Ms series 




k = 1 


0.0381237 


0.04(1 + 0.05) 


0.04017(205] 


k = 2 


0.0382058 


0.04(1 + 0.05 - 0.02) 


0.03955[71] 


k = 3 


0.0382161 


0.04(1 + 0.054 - 0.019 - 0.004) 


0.03939[17] 


Order 


a 


TZ series 


n 


k = 1 


0.0414570 


0.04(1 - 0.02) 


0.04043[103] 


k = 2 


0.0394420 


0.04(1-0.01 + 0.01) 


0.03944(47] 


k = 3 


0.0391507 


0.04(1 + 0.003 + 0.002 + 0.001) 


0.03941 [4] 



Table 1. Results for 7^, the QCD corrections to Re+e~, in {k + I)*'' order (N'^LO) at 
an energy Q/A^s = 340. The upper and lower sub-tables list, respectively, the MS 
and optimized results. The columns give the couplant value, the rough form of the 
series, and the result for TZ with an error estimate corresponding to | VkO,^^^ |, the 
magnitude of the last term included in the perturbation series. 



0.043 r 1 

0.042^ [ z 

0.041 \ I { 

0.039^ n *' ^ 

0.038^ ^ ^ 

0.037 E 

1 2 3 

k 

Figure 1: Results for TZ, the QCD corrections to R^+^-, in (k + 1)*'^ order (N'^LO) at an energy 
Q/Ams = 340. The larger points displaced leftwards are in the MS scheme, while the smaller 
points displaced rightwards are the optimized results. In both cases the error bars correspond 
to |rfca''+^ |. 
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Example 2: Q/Ams = 68. 



Order 




T^Ms series 




k = 1 


0.0507097 


0.05(1 + 0.07) 


0.05433[362] 


k = 2 


0.0509032 


0.05(1 + 0.07-0.03) 


0.05287(169] 


k = 3 


0.0509356 


0.05(1 + 0.07- 0.03 - 0.01) 


0.05236[54] 


Order 


a 


TZ series 


n 


k = 1 


0.0568587 


0.06(1 - 0.03) 


0.05496[190] 


k = 2 


0.0525541 


0.05(1 - 0.02 + 0.02) 


0.05256[112] 


k = 3 


0.0520416 


0.05(1 + 0.002 + 0.003 + 0.002) 


0.05245[13] 



Table 2. Results (as Table 1) for Q/Aws = 68. 



0.058 

0.056 ; 

0.054 ; 

0.052; 

0.050- 

1 2 3 

k 

Figure 2: As Fig. 1 but for Q/Ay^ = 68. 
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6.3 Low-energy examples 

Next we turn to lower-energy examples, where the differences between MS and OPT are more 
dramatic. With nj = 2 the f3 function's leading coefficients are 

b=-, c=—, (6.13) 

and da EI] 

48241 

= ^ = 5,77598, (6.14) 

18799309 68881 

^ 902016 12528^ ^ ^ ^ 

The TZ coefficients, in the MS scheme, are [IHl EOj [6] 

= 1.755117, rf^ = -9.14055, rf^ = -123.18799. (6.16) 

(Again, the exact values [6j were used in our calculations.) Inserting these values in Eq. ()2.40p 
yields 

P2 = -9.92498, p3 = -115.21021. (6.17) 

Tables 3-7 and Figures 3-7 give results at successively lower energies; Q/Ams = 5, 2, 1.7, 1.5 
and 0. One sees in the MS results the characteristic symptoms of an asymptotic series; after 
initially seeming to converge, the series starts to go bad, with the error estimate increasing with 
order. In Example 3 the effect is just visible in the k = 3 result, but in Examples 4 and 5 
the effect becomes more dramatic. In Example 6 there is no /c = 3 MS result at all since for 
Q/A^g < 1.645 there is no solution to the k = 3 int-/3 equation. At still lower values of Q/A^g, 
below 1.396 and 1 respectively, the k = 2 and k = 1 MS int-/3 equations have no solution. 

By contrast, the optimized results show a monotonic decrease in the expected error at higher 
orders. The k = 1 results, in Examples 5 and 6 particularly, are very uncertain at low energies — 
indeed, for Q/A-^ < 1.438 there is no solution to the k = 1 optimal int-/3 equation, Eq. (I3.14p . 
However, for k = 2 and 3 the optimized results improve very significantly. Indeed, for k = 2 
and 3 there are optimized results down to zero energy, because the optimized (3 function turns 
out to have a non-trivial fixed point (see Appendix B). The new A; = 3 results provide solid 
confirmation of the earlier k = 2 results [21 HI [5] . As shown in the table for example 7, the Q = 
result for TZ improves from 0.3 it 0.3 to 0.2 it 0.1. 
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Example 3: Q/Ams = 5. 



Order 




T^Ms series 




k = 1 


0.0862557 


0.09(1 + 0.15) 


0.099[13] 


k = 2 


0.0902494 


0.09(1 + 0.16 -0.07) 


0.098[7] 


k = 3 


0.0911287 


0.09(1 + 0.16 - 0.08 - 0.09) 


0.090[8] 


Order 


a 


TZ series 


n 


k = 1 


0.117285 


0.12(1 - 0.09) 


0.106[11] 


k = 2 


0.0952429 


0.10(1 - 0.05 + 0.05) 


0.095 [5] 


k = 3 


0.0899359 


0.09(1 - 0.01 - 0.02 + 0.04) 


0.091 [4] 



Table 3. Results for Q/A^g = 5. 



0.12; 

0.11; I 

0.10 O T 

(I 

0.09; ^ 

0.08 E 

1 2 3 
k 

Figure 3: As Fig. 1 but for Q/Ams = 5. 
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Example 4: Q/A^g = 2. 



Order 




T^Ms series 




k = 1 


0.1626471 


0.16(1 + 0.29) 


0.209[46] 


k = 2 


0.1963533 


0.20(1 + 0.34 -0.35) 


0.195 [69] 


k = 3 


0.2193679 


0.22(1 + 0.39 -0.44 - 1.30) 


-0.08 + 0.29 


Order 


a 


1Z series 


H 


k = 1 


0.3648099 


0.36(1-0.21) 


0.288[77] 


k = 2 


0.1725913 


0.17(1-0.17 + 0.18) 


0.173[31] 


k = 3 


0.1421756 


0.14(1 -0.08 -0.09 + 0.17) 


0.143[25] 



Table 4. Results for Q/Ams = 2. 



1 2 3 

k 

Figure 4: As Fig. 1 but for Q/A^ = 2. (The k = 3 MS result is slightly negative; only its error 
bar is visible.) 



0.3 
0.2 
0.1 
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Example 5: Q/A^ = 1.7. 



Order 




T^Ms series 




k = 1 


0.1966624 


0.20(1 + 0.35) 


0.265 [68] 


k = 2 


0.2691684 


0.27(1 + 0.47-0.66) 


0.218[178] 


k = 3 


0.4153849 


0.42(1 + 0.73 - 1.58 - 8.83) 


-3.60 + 3.67 


Order 


a 


1Z series 


H 


k = 1 


0.6669931 


0.67(1 -0.28) 


0.477[190] 


k = 2 


0.1970393 


0.20(1 -0.24 + 0.25) 


0.199[49] 


k = 3 


0.1530735 


0.15(1-0.11-0.11 + 0.22) 


0.153[34] 



Table 5. Results for Q/Ams = 1-7. 



I : 

1 2 3 

k 

Figure 5: As Fig. 1 but for Q/A^g = 1.7. (Only the tip of the huge error bar for the k = 3 MS 
result is visible.) 



0.6 
0.4 
0.2 
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Example 6: Q/A-^ = 1.5. 



Order 




T^Ms series 




k = 1 


0.236877 


0.24(1 + 0.42) 


0.335[98] 


k = 2 


0.431322 


0.43(1 + 0.76 - 1.70) 


0.02 + 0.73 


k = 3 


no solution 






Order 


a 


1Z series 


n 


k = 1 


2.4690661 


2.5(1 - 0.42) 


1.4 + 1.0 


k = 2 


0.2173977 


0.22(1 -0.31 + 0.33) 


0.221[71] 


k = 3 


0.1605183 


0.16(1-0.13-0.13 + 0.26) 


0.161[42] 



Table 6. Results for Q/Ams = 1-5. 



1.5 



1.0 



0.5 



0.0^ 



Figure 6: As Fig. 1 but for Q/A^s = 1-5. (There is no = 3 MS result in this case. 
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Example 7: Q/A^ = (fixed point). 



Order 


a 


IZ series 


n 


k = 1 


no solution 






k = 2 


0.2635259 


0.26(1 -0.76 + 1.01) 


0.330[267] 


k = 3 


0.1800794 


0.18(1 -0.25 -0.16 + 0.44) 


0.185[79] 



Table 7. Results for the infrared fixed-point limit, Q/A^g = 0. There are no MS 
results in this case. 



0.8 [ 

0.6 

0.4 

0.2 

0.0 



Figure 7: As Fig. 1 but for the infrared fixed-point limit, Q/Ams = 0. There are no MS results 
in this case, and no /c = 1 optimized result. 
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7 Concluding Remarks 



We first summarize the lessons of tlie numerical examples, which compared the MS(;U = Q) 
and optimized results. At moderately high energies, the differences are well within the error 
estimates; the main advantage of optimization here is to achieve greater precision. At low 
energies, however, the optimized results show steady convergence, while the MS results begin 
to show the typical pathologies of a non-convergent asymptotic series. We can expect these 
pathologies to eventually show up in the MS results at higher energies when the series is taken 
to high enough order. In this sense the low-energy examples are a "preview" of the divergent- 
series problems to be expected in MS or any fixed RS. 

Whether or not the optimized results ultimately converge is a matter of conjecture at present. 
However, the lessons of toy models [221 [23] and from the linear 5-expansion for the anharmonic 
oscillator [15], together with the present results, suggest that it is a very real possibility. We 
should say at once that we would not expect the optimized results to converge to the exact, 
physical answer. We certainly expect there to be nonperturbative contributions (higher- twist 
terms that involve powers of e~^/ ^^"^ and hence are invisible to perturbation theory) . A physically 
defined quantity TZ in QCD is, in general, a sum of perturbative and nonperturbative terms, 
T^cTt + '^nonpert- The issue, though, is whether this decomposition can be made unique and 
physically meaningful or whether it is inherently ambiguous and dependent on RS. It seems 
reasonable that T^nonpcrt should involve only infrared physics and so should be calculable, in 
principle, without any essential need for renormalization. If that is so, there would have to be 
a version of perturbation theory that converges to a RS-invariant result, T^pert- Optimization, 
where at each order we demand that the result be "as RS invariant as possible," is the natural 
candidate to be that version of perturbation theory — without any need to invoke extra tricks 
(such as Fade approximants, Borel summation, etc.). 

The good convergence of the optimized results in the e^e^ case remains true even down to 
zero energy, where the third-order finding [3ll31[5] of a limit TZ = 0.3 it 0.3 is confirmed and made 
more precise; TZ = 0.2 it 0.1. This result is very important because there are many indications 
from phenomenology (see the mini-review in Ref. [5]) that the QCD couplant does freeze at low 
energies. Usually freezing is something that theorists put in by hand, but here it is an outcome, 
a prediction. There is nothing in the optimization approach that forces freezing to occur; the 
fact that it does for R^+e- is due to the P2,P3 values resulting from the Feynman-diagram 
calculations. 

The methods described here could be applied to many other perturbative physical quantities. 
In particular, fourth-order results are now available for r-lepton and decay widths and for 
some deep-inelastic-scattering sum rules [H [36]. Frevious studies of OFT applied to these 
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quantities [31 [371 [3E] could now be extended to fourth order. 

We have not discussed here physical quantities that explicitly involve parton distribution 
functions or fragmentation functions. Such quantities are plagued by another kind of ambi- 
guity; factorization-scheme dependence — an ambiguity similar to, and entangled with, RS 
dependence. The "principle of minimal sensitivity" can be applied here, too. Unfortunately, 
the original analysis [39] was stymied by an algebraic error, corrected later in Ref. [IQ] (see also 
[41j). These papers use the language of structure- function moments — which, while natural the- 
oretically, is perhaps not very convenient phenomenologically. A purely numerical approach to 
"optimization" [421 US] is certainly feasible, but laborious. It would be valuable to somehow re- 
formulate Ref. [40] 's results in a way that would combine easily with practical methods for using 
and empirically determining parton distribution functions. We suspect, based on the important 
results of |421 143j. that several important QCD cross sections are currently underestimated, and 
that "optimization" could significantly reduce the theoretical uncertainties of many others. 
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Appendix A: Proof of the "complete sum" identity 

This appendix provides a proof of the identity, mentioned in subsection 14.31 for the "complete 



^ cja^ B,+,{a) = 1 i = {l),2,...,k. (A.l) 

From the definition of the Pj{a) and Bj{a) functions we have Eqs. (j2.19p and (j2.20p . Sub- 
stituting into the l.h.s. of (|A.ip yields 



i=o 







B{xY ■ 

Pulling the integration outside the sum and re-grouping leads to 



The integrand can now be recognized as 



,^T,r ^ dB{x)\ d ( x'~^ . ,^ ^. 



B{x)^ \ dx J dx \B{x) J 
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so the integral can be done immediately, giving 



Bja) 

ni-l 



B{x) 



(A.5) 



x=0 



For i > 1 the lower endpoint (x = 0) makes no contribution, so the result is unity, as claimed. 

The case i = 1 requires special treatment. As in the definition of Hi{a), Eqs. (j4.3l I4.4p . we 
should interpret the i = 1 case of Eq. (jA.ip above as 

k 

B{a) - ^Cja^Bj+i{a) = 1. (A.6) 
i=i 

The proof parallels the steps above except that there is now a contribution from the x = 
endpoint that results in the cancellation of the B{a) term. 



Appendix B: Infrared limit and fixed-point behaviour 

In asymptotically free theories, perturbation theory works best at high energies. By investigating 
low energies we can learn important lessons about how perturbation theory can go bad, and why 
RS choice is crucial. As the physical scale Q is lowered, and the effective couplant grows, we can 
expect the physical quantity IZ to either (z) go to infinity at some finite energy of order A, or 
(ii) tend to a finite limit as Q — )■ 0. The latter scenario is usually said to happen if and only if 
the /? function has a zero at some finite a = a*, called a "fixed point." That statement, however, 
is too naive because "the /? function" is not a unique object; it depends on RS. 

At second order in QCD, for < 8, the /3 function has no non-zero fixed point in any RS. 
At higher orders, though, the question depends entirely on the RS choice. In the MS scheme the 
C2 and C3 coefficients are positive, so no fixed point exists and MS results, at third and fourth 
orders, go to infinity at some Q of order Ams- However, the optimization procedure does give 
finite results at arbitrarily low energies in the e'^e~ case. At third order [5] TZ{e'^e~) approaches 
a finite limit, O.SitO.S. At fourth order we find 0.2ib0.1. These results can be found by applying 
the algorithm of Section [5] at lower and lower Q, but can also be obtained much more simply 
because [24] the optimization equations greatly simplify at a fixed point. The results of Ref. |24j 
(converting notation: p'^'^ = p2 — \(? and p'^'^ = \p^) are as follows. At third order the optimized 
a* is given by the smallest root of the quadratic equation 

-^+ca* + 2,{p2--/^a*^ = (B.l) 
and the limiting value of 'RP'^ is 
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At fourth order the corresponding equations are 

§ + T?^"* + 7'52a*' + 2p3a*' = 0, (B.3) 
d4 lb 4 

and 

-r>*(4) * ( 249 13 ^ 1 , /T3 ^\ 

^ V^ + 64^" +16^^" J- ^^-'^ 

The fixed-point limit of Eq. (j4.13p . the new formula for the Sm = ( "^p^ ) coefficients, can 
be written as 

Sm = _ [{k -2m- l)tm -{k- 2m)im-i\ , (B.5) 

or, equivalently, 

•Sm = _ [{k - 2m)cm - im] , (B.6) 

where Sm = SmCL*"^, and dm = CmCL*"^, and tm is a partial sum of /3-function terms: 

m 

^m = 5^cX''. (B.7) 

i=0 

These new formulas simplify the task of generalizing Ref. |24J's results to higher orders. 

The occurrence of a fixed point is not inevitable in OPT; the a* equation, (jB.ip or ()B.3p . 
may or may not have a positive, real root. A small positive a* is found for TZ{e~^e~), for nj = 2, 
because the invariants p2 and ps are negative and sizeable. Interestingly, for nj above about 6 
(depending on the assumed electric charges of the extra quarks) one does not find a solution to 
Eq. (jB.3p . However, a finite infrared limit still exists, but it occurs by a new mechanism; we 
hope to report on this phenomenon in a future publication. 
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